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The simplest unified extension of the Minimal Supersymmetric Standard Model 
with bilinear R-Parity violation provides a predictive scheme for neutrino masses 
which can account for the observed atmospheric and solar neutrino anomalies. De¬ 
spite the smallness of neutrino masses R-parity violation is observable at present 
and future high-energy colliders, providing an unambiguous cross-check of the 
model. 

1. Introduction 

The announcement of high statistics atmospheric neutrino data by the 
SuperKamiokande collaboration 1 has confirmed the deficit of muon neutri¬ 
nos, especially at small zenith angles, opening a new era in neutrino physics. 
Although there may be alternative solutions of the atmospheric neutrino 
anomaly 2 it is fair to say that the simplest interpretation of the data is in 
terms of v tJ to v T flavor oscillations with maximal mixing. This excludes a 
large mixing among v T and v e , in agreement also with the CHOOZ reactor 
data 3 . On the other hand the persistent disagreement between solar neu¬ 
trino data and theoretical expectations 4 has been a long-standing problem 
in physics. Recent solar neutrino data 5 are consistent with both vacuum 
oscillations and MSW conversions. In the latter case the large mixing angle 
solutions are now clearly preferred 2 . 

Many attempts have appeared in the literature to explain the data. Here 
we review recent results 6 obtained in a model 7 which is a simple extension 
of the MSSM with bilinear R-parity violation (BRpV). This model, despite 
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being a minimal extension of the MSSM, can explain the solar and atmo¬ 
spheric neutrino data. Its most attractive feature is that it gives definite 
predictions for accelerator physics for the same range of parameters that 
explain the neutrino data. 

2. Bilinear R-Parity Violation 

2.1. The Model 

Since BRpV SUSY has been discussed in the literature several times 6, ' ,s 
we will repeat only the main features of the model here. We will follow the 
notation of 6 . The simplest bilinear Jji v model (we call it the MSSM) is 
characterized by three additional terms in the superpotential 

W = Wmssm + W^ p (1) 

where Wmssm is the ordinary superpotential of the MSSM and 

W$ p = eiUHu. ( 2 ) 

These bilinear terms, together with the corresponding terms in the soft 
SUSY breaking part of the Lagrangian, 

Csoft = + BietLiHu (3) 

define the minimal model, which we will adopt throughout this paper. The 
appearance of the lepton number violating terms in Eq. (3) leads in general 
to non-zero vacuum expectation values for the scalar neutrinos (£*), called 
Vi in the rest of this paper, in addition to the VEVs vjj and Vn of the 
MSSM Higgs fields H ° and H Together with the bilinear parameters 
£i the Vi induce mixing between various particles which in the MSSM are 
distinguished (only) by lepton number (or R-parity). Mixing between the 
neutrinos and the neutralinos of the MSSM generates a non-zero mass for 
one specific linear superposition of the three neutrino flavor states of the 
model at tree-level while 1-loop corrections provide mass for the remaining 
two neutrino states, see 6 . 

3. Neutrino Masses and Mixings 

3.1. Tree Level Neutral Fermion Mass Matrix 

In the basis i/j ot = (—iA', —?'A 3 , R], F e , ^ r ) the neutral fermions 

mass terms in the Lagrangian are given by 

C m = M N ip° + h.c. 


( 4 ) 
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where the neutralino/neutrino mass matrix is 

M n = 
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where a,i = (— ^g'vii \q v u 0, £*). This neutralino/neutrino mass matrix is 
diagonalized by 

AT*M N A 1 = diag(m % o, m x o,m x o,m x o,rn vl ,m V2 ,m V3 ) 


(7) 


3.2. Approximate Diagonalization at Tree Level 


If the IjLp parameters are small, then we can block-diagonalize Mn approx¬ 
imately to the form diag(ra e //, AA x o) 


m e ff = -m ■ Mj;m T 


Mig 2 + M 2 g ' 2 
Adet ( A 4 x o ) 


1 A 2 AeA^AeA^ 
A e A M A“ A^A-t- 
y A e A r A m A t A t J 


( 8 ) 


The matrices N and V v diagonalize A4 x o and m e // 

N*A4 x oN'=diag(m x 9 ) ; Vj m e f fV v = diag (0,0, m„), (9) 


where 


m v = Tr(m e ff) 


Miff 2 + Mig' 2 ^ 2 

4det(Ai x o) 


( 10 ) 


3.3. Approximate Formulas for 1 Loop 

3.3.1. The masses 


Looking at the numerical results 6 we found that the most important con¬ 
tribution came from the bottom-sbottom loop. To gain an analytical un¬ 
derstanding of the results we expanded the exact results in the small ftp 
parameters. The result is 9 


M v ~ cq 
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( 11 ) 
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where 


M l9 2 + M 2 g' 2 
4 clet(M x o) 


ci = 
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Diagonalization of the mass matrix gives 9 


( 12 ) 
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(13) 

(14) 

(15) 


The formula for m „ 3 is the tree-level formula that we used to hx the scale 
of the atmospheric neutrinos by choosing |A|. Details of the derivation can 
be found in Ref. 9 where the second most important contribution, coming 
from the loop with charged Higgs/charged leptons, is also discussed. 


3.3.2. The mixings 

The atmospheric angle is easily obtained in terms of the ratio A 2 /A 3 . For 
the solar angle in the same approximation we also get a simple formula 9 , 

2e 2 


tan 0sol — / 

(£2 + e 3 ) 2 

that is also in very good agreement with the exact result. 


(16) 


4. Results for the Solar and Atmospheric Neutrinos 
4.1. The masses 

The BRpV model produces a hierarchical mass spectrum for almost all 
choices of parameters. The largest mass can be estimated by the tree level 
value using Eq. (15). Correct A m 2 trn can be easily obtained by an appro¬ 
priate choice of |A|. The mass scale for the solar neutrinos is generated 
at 1 -loop level and therefore depends in a complicated way in the model 
parameters. However, in most cases the result of Eq. (14) is a good ap¬ 
proximation and there is no problem in having both A m 2 tm and A m 2 olar 
set to the correct scales. 


4.2. The mixings 

Now we turn to the discussion of the mixing angles. We have found that if 
e 2 /|A| -C 100 then the 1-loop corrections are not larger than the tree level 
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results and the flavor composition of the 3rd mass eigenstate is approxi¬ 
mately given by 

U a 3 « A a /|A| (17) 

As the atmospheric and reactor neutrino data tell us that —•> v T oscilla¬ 

tions are preferred over v tJ —* v e) we conclude that 

A e < A m ~ A t (18) 

are required for BRpV to fit the data. This is sown in Fig. 1 a). We cannot 
get so easily maximal mixing for solar neutrinos, because in this case U e 3 
would be too large contradicting the CHOOZ result as shown in Fig. 1 b). 



0.5 1 

ivi/yA^ 


0.5 1 

Figure 1 . a) Atmospheric angle as a function of |A^|/vAj + A?. b) U g 3 as a 


function of |A e |/y Aj) + A?. 


We have then two scenarios. In the first one, that we call the mSUGRA 
case, we have universal boundary conditions of the soft SUSY breaking 
terms. In this case we can show 6 that 


Ae 

A„ 


(19) 


Then from Fig. 1 b) and the CHOOZ constraint on t/g 3 , we conclude that 
both ratios in Eq. (19) have to be small. Then from Fig. 2 we conclude that 
the only possibility is the small angle mixing solution for the solar neutrino 
problem. In the second scenario, which we call the MSSM case, we consider 
non-universal boundary conditions of the soft SUSY breaking terms. We 
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have shown that even a very small deviation from universality of the soft 
parameters at the GUT scale relaxes this constraint. In this case 



( 20 ) 


Then we can have at the same time small t/ 2 3 determined by A e /A M as in 
Fig. 1 b) and large tan 2 (0 so j ar ) determined by Ce/e^ as in Fig. 2 b). 



5. Probing Neutrino Mixing via Neutralino Decays 

If R-parity is broken, the neutralino is unstable and it will decay through 
the following channels: Xi v i v 3 v ki qq , Vi lj~ l ifqq', Vi 7 . It was 
shown a in Ref. 10 , that the neutralino decays well inside the detectors and 
that the visible decay channels are quite large. This was fully discussed 
in Ref. 10 where it was shown that the ratios | A.^ /Aj | and | /ey | were 
very important in the choice of solutions for the neutrino mixing angles. 
What is exciting now, is that these ratios can be measured in accelerator 
experiments. In the left panel of Fig. 3 we show the ratio of branching ratios 
for semileptonic LSP decays into muons and taus: BR{\ —► Hq'q)/BR(x —► 
rq'q) as function of tan 2 d atm . We can see that there is a strong correlation. 
The spread in this figure can in fact be explained by the fact that we do not 
know the SUSY parameters. This is illustrated in the right panel where we 


a The relation of the neutrino parameters to the decays of the neutralino has also been 
considered in Ref. 11 . 
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considered that SUSY was already discovered with the following values for 
the parameters, 


M 2 = 120 GeV, /i = 500 GeV ,, tan/3 = 5, m 0 = 500 GeV. A = -500 GeV (21) 


ed 

PQ 


g °- 5 


0.1 




Figure 3. Ratios of semileptonic branching ratios as functions of tan Q a tm- On 
the left for random SUSY values and on the right for the SUSY point of Eq. (21) 


6. Probing Neutrino Mixing via Charged Lepton Decays 

If R-parity is broken the lightest supersymmetric particle (LSP) will decay. 
If the LSP decays then cosmological and astrophysical constraints on its 
nature no longer apply. Thus, within R-parity violating SUSY a priori any 
superparticle could be the LSP. We have studied 12 the case where a charged 
scalar lepton, most probably the scalar tau, is the LSP. The production and 
decays of f, as well as the decays of e and p,, and demonstrate that also 
for the case of charged sleptons as LSPs neutrino physics leads to definite 
predictions of various decay properties. This is shown in Figs. 4 and 5 

7. Conclusions 

The Bilinear R-Parity Violation Model is a simple extension of the MSSM 
that leads to a very rich phenomenology. We have calculated the one-loop 
corrected masses and mixings for the neutrinos in a completely consistent 
way, including the RG equations and correctly minimizing the potential. 
We have shown that it is possible to get easily maximal mixing for the 
atmospheric neutrinos and both small and large angle MSW. 
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Figure 4. a)e + e —> ll production croi 
0.8TeV, b) Charged slepton decay length 



100 150 200 250 300 350 400 


m ii t GeV ] 

section at a Linear Collider y/s = 
a linear collider with yfs = 0.8TeV. 



Figure 5. Ratios of branching ratios for scalar tau decays versus tan 2 #Q. The 
left panel shown all data points, the right one refers only to data points with 
£ 2/^3 restricted to the range [0.9,1.!]. 


We emphasize that the LSP decays inside the detectors, thus leading to 
a very different phenomenology than the MSSM. The LSP can be either the 
lightest neutralino, like in the MSSM, or a charged particle, must probably 
the lightest stau. In both cases we have shown that ratios of the branching 
ratios of the LSP can be correlated with the neutrino parameters. 

If the model is to explain solar and atmospheric neutrino problems many 
signals will arise at future colliders. These will probe the neutrino mixing 
parameters. Thus the model is easily falsifiable! 
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